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Abstract: Using gauge-gravity duality, we extend thermodynamic studies and present 
results for thermal screening masses in strongly coupled J\f = 2* supersymmetric Yang- Mills 
theory. This non-conformal theory is a mass deformation of maximally supersymmetric = 4 
gauge theory. Results are obtained for the entropy density, pressure, specific heat, equation 
of state, and screening masses, down to previously unexplored low temperatures. The tem- 
perature dependence of screening masses in various symmetry channels, which characterize 
the longest length scales over which thermal fluctuations in the non-Abelian plasma are cor- 
related, is examined and found to be asymptotically linear in the low temperature regime. 
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1 Introduction 

A fundamental characteristic of any system in thermal equilibrium is the correlation length: 
the longest distance scale over which spatial fluctuations are significantly correlated. More 
precisely, the correlation length describes the slowest exponential fall-off of correlation func- 
tions at asymptotically large separation. One may choose to restrict attention to correlators 
of operators with specified symmetries, and define correlation lengths in individual symmetry 
channels. Inverse correlation lengths have units of energy and are commonly referred to as 
screening masses. These describe the long distance fall-off of correlations in fluctuations with 
specified quantum numbers, and hence characterize the nature of infrared effective degrees of 
freedom. 

In, for example, a relativistic QED plasma, Debye screening leads to exponential fall-off 
of the electric field induced by an external test charge, (-E(x)) oc e~l^l/^, and corresponding 
behavior in the correlator of electric field fiuctuations. The inverse correlation length, or 
Debye screening mass, mr, = 1/^, is 0{eT). Fermionic fiuctuations in a weakly coupled 
relativistic plasma have correlations which decrease exponentially with distance on a shorter 
0(1/T) length scale, corresponding to the inverse of the lowest fermionic Matsubara frequency. 
Static magnetic fields are not screened at long distances, so the associated correlation length 
is infinite. 

In a non-Abelian QCD plasma, one may isolate the physics of color-electric screening in 
a gauge-invariant and non-perturbative manner by examining correlators of operators which 
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are odd under Euclidean time reflection [1]. The Debye screening length can be defined as the 
longest correlation length in time-reflection odd symmetry channels.^ At asymptotically high 
temperatures the Debye screening length may be calculated perturbatively. For an SU{Nc) 
gauge group with A'^f Dirac fermions one finds, at leading order, ttid = [^^c + gA^f] ^^"^ g{T)T. 
Long distance properties of QCD at high temperature are effectively described by three- 
dimensional pure Yang-Mills gauge theory with a dimensionful coupling = g{T)'^T. This 
theory exhibits three-dimensional confinement and the generation of a non-perturbative 0(53) 
mass gap. Therefore, the longest distance correlations in asymptotically high temperature 
QCD are associated with static magnetic fiuctuations, and have an 0[{g(T)'^T)^^] correlation 
length. 

At non-asymptotic temperatures, reliable calculations of QCD screening masses require 
numerical lattice simulations.^ At temperatures of a few times the deconfinement tempera- 
ture, T/Tc ~ 1-4, which is the temperature range probed by heavy ion collisions, the dynamics 
of QCD plasma is strongly coupled. Evidence for the strongly coupled nature of QCD plasma 
in this regime includes the results for screening masses obtained from lattice simulations (see 
ref. [3] for a review) , as well as the success of low viscosity hydrodynamic simulations in repro- 
ducing the collective fiow observed at RHIC, and the observed suppression of high transverse 
momentum jets in RHIC collisions [4]. 

Numerical lattice simulations are limited to equilibrium Euclidean observables; with very 
limited exceptions, non-equilibrium Minkowski space dynamics cannot be extracted reliably 
from lattice simulations. Consequently, there has been much interest in studying dynamic 
properties of non-Abelian gauge theories which mimic some aspects of QCD and to which the 
techniques of gauge/gravity (or AdS/CFT) duality [5-7] may be applied. (See, for example, 
refs. [8-16].) Gauge/gravity duality reformulates the strong coupling, large Ac dynamics 
of a suitable gauge theory in terms of classical supergravity in an asymptotically anti-de 
Sitter spacetime. The simplest, and best understood, example is maximally supersymmetric 
SU{Nc) Yang-Mills theory (M = 4 SYM). Although this theory contains adjoint representation 
fermions and scalars not present in QCD, at non-zero temperature the resulting non-Abelian 
plasma shares many qualitative similarities with the quark-gluon plasma of QCD. Screening 
masses obtained via gauge-gravity duality for strongly coupled (A = g'^Nc ^ 1) = 4 SYM, in 
the A^c — ^ 00 limit, were discussed in refs. [17, 18]. The conformal invariance of AA = 4 SYM 
implies that thermal screening masses are exactly proportional to T. Ref. [17] compared 
the values of screening masses of strongly coupled AA = 4 SYM with results of lattice QCD 
simulations at T ~ 2Tc. In QCD, this temperature lies within the window 1.5 Tc < T < 4Tc 
where screening masses are well described as scaling linearly with T [19, 20]. Ratios of 
screening masses in different symmetry channels compare rather well between QCD and AA = 4 

^Euclidean time reflection Rr corresponds to the product TC of time reversal and charge conjugation. 
This definition of the Debye screening length applies only to TC-invariant equilibrium states in TC-invariant 
theories. 

^For example, up to temperatures as high as 10^ Tc, the non-perturbative O(g^T) correction to the Debye 
screening mass in QCD is larger than the leading 0{gT) perturbative contribution [2]. 
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SYM but, in absolute size, the J\f = 4: SYM screening masses (divided by T) are roughly a 
factor of two larger than in QCD [17]. 

It is not currently known whether quantitative differences between QCD and = 4 SYM 
plasma properties are dominantly due to (i) the differing constituents of the non-Abelian 
plasmas, [ii) the comparison of QCD at T /Tq ~ 1.5-4, where it is neither weakly coupled nor 
infinitely strongly coupled, with the A ^ 1 limit of SYM instead of some 0(1) value of A, or 
{in) the comparison with = oo SYM instead oi = 3. Available lattice results suggest 
that the dependence of plasma thermodynamics and screening masses is quite mild [19, 21], 
so possibility [in) is unlikely to be dominant. Possibility {ii) could be studied by computing 
1/A corrections to the leading strong coupling behavior described by classical supergravity. 
However, for many observables this requires knowledge of higher derivative corrections to the 
supergravity action which have not yet been fully worked out. The possibility most accessible 
to study is (z): the dependence of various observables on the particular constituents of a 
non-Abelian plasma. Part of the reason that screening masses are larger in A/" = 4 SYM than 
in QCD is certainly the fact that there are more active degrees of freedom contributing to 
screening in an AA = 4 SYM plasma — adjoint representation fermions and scalars instead 
of fundamental representation quarks. One would expect that reducing the number of light 
degrees of freedom would yield screening masses which are closer to those of a QCD plasma. 
This may be explored, in strongly coupled plasmas, by studying theories with known gravity 
duals which are more QCD-like than AA = 4 SYM. Examining non-conformal deformations of 
AA = 4 SYM will shed light on the sensitivity of ratios of screening masses to departures from 
conformal invar iance. 

In this paper, we explore these issues by computing screening masses and thermodynam- 
ics of a mass-deformation of maximally supersymmetric Yang- Mills theory known as AA = 2* 
SYM [22-24]. This theory differs from A/" = 4 SYM by the addition of mass terms and in- 
teractions for an adjoint representation hypermultiplet, in a manner which preserves N = 2 
supersymmetry. This theory has one new free parameter, the deformation mass m. The 
theory is non-conformal, with a nontrivial renormalization group fixed point (with adjustable 
coupling A) at asymptotically high energies. An intrinsic strong scale A appears, analogous 
to AqcD) but for large A this scale is comparable to the mass deformation, A ~ m [22]. 

The study of M = 2* gauge theory at non-zero temperature using holographic methods 
was pioneered by Buchel, Liu, and collaborators [24-26]. Various properties of AA = 2* plasma 
have been studied, including shear and bulk viscosities [25, 27, 28], as well as heavy quark 
drag, momentum broadening, and jet quenching [29]. Reference [26] laid the groundwork for 
numerically constructing the supergravity background dual to the thermal equilibrium state 
of AA = 2* SYM, and mapped out the dictionary that translates between supergravity fields 
and observables in the quantum field theory. Of direct relevance to our work is their result for 
the free energy density in the range < m/T < 6.'^ Noticeable deviation from the behavior 
of AA = 4 plasma was estimated to occur near the upper end of this range. 

^In the later work [27], hydrodynamic coefficients were evaluated up to m/T ~ 12. 
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We extend the results of ref. [26] on = 2* thermodynamics to substantially lower tem- 
peratures, m/T < 33, and provide evidence that this data allows one to probe the asymptot- 
ically low temperature regime. We argue that, despite the presence of a non-zero mass scale 
m, this theory (in the large and large A limit) has no thermal phase transition at any 
non-zero temperature. 

Our primary goal is the evaluation of thermal screening masses for multiple symmetry 
channels throughout the temperature range < m/T < 33. We extrapolate the screening 
mass curves to vanishingly small temperatures, T/m — t- 0, where components of the heavy 
hypermultiplet decouple from the dynamics due to exponential Boltzmann suppression, and 
extract results for the limiting plasma which contains only light degrees of freedom. 

The paper is structured as follows. In Sec. 2, we briefly review the field theoretic for- 
mulation of A/" = 2* gauge theory and its gravitational representation. In Sec. 3, we discuss 
equilibrium thermodynamic quantities such as the entropy density, pressure, specific heat, and 
equation of state. We compute screening masses from linearized fluctuations of supergravity 
modes in Sec. 4. In Sec. 5, we compare the screening masses of = 2* SYM to those of A/" = 4 
SYM and to available lattice QCD data. The holographic renormalization and derivation of 
thermodynamic quantities from the gravitational dual are discussed, at some length, in Ap- 
pendix A. The numerical method we used to find the background geometry is explained in 
Appendix B, and fluctuations of the gravitational fields we use to compute screening masses 
are described in Appendix C. Additional material related to the appendices can be found 
online [30]. 

2 J\f = 2* gauge theory and its gravity dual 

M = 2* supersymmetric Yang-Mills theory is a mass deformation of A^ = 4 SYM. We consider 
the theory on M}'^ with gauge group SU (Nc) in the 't Hooft limit: Nc ^ oo and 5^ — >• with 
A = g'^Nc held fixed. If the A/" = 4 field content is grouped into A/" = 1 superfields, there is a 
vector multiplet and three adjoint chiral multiplets {i = 1,2,3). The Lagrange density of 
M = 2* SYM is obtained by adding the term 



to the A/" = 4 SYM Lagrange density Cj\f = ii. Appropriate field redefinitions can remove any 
phase in the mass, so m may be taken real and positive without loss of generality. The two 
chiral fields (<I>i, <I>2) comprise a single AA= 2 massive hypermultiplet, and the remaining fields 
form an AA = 2 vector multiplet. Once auxiliary field constraints have been solved, the above 
contribution to the superpotential induces conventional mass terms for two Weyl fermions 
and two complex scalars, as well as terms trilinear in the scalars. In order to identify the 
gravity dual of the M = 2* SYM theory, it is convenient to write the mass deformation in 
terms of relevant operators in irreducible representations of the AA = 4 i?-symmetry group 
SO{Q)ii. We can distinguish two contributions, first there is an irreducible scalar mass 




(2.1) 




(2.2) 
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This operator, which is an element of the 20' of S0{6), appears in the 4D (Euchdean) J\f = 2* 
action Sj\f = 2* as —{2/g^) f d^xm?02- Hence, it contributes the usual positive bosonic mass 
terms for the hypermultiplet scalars and 02- But it also destabilizes the scalar 03 belonging 
to the vector multiplet. The second contribution introduces a mass for the fermions in the 
hypermultiplet 

O3 = tr - ^H'Pu 4>\] - ^H'PI^ H + h.c.) - I mtr {\(t>i\'' + |02|' + , (2-3) 

and appears in Sv = 2* as —{2/g^) f df^xmO^. One may obtain the expression for by first 
deriving the on-shell Lagrange density for mass-deformed = 4 SYM in terms of component 
fields, isolating the m-dependent terms, and subtracting m? O2 from them. Note that is 
modified, from its massless limit, by the £'0(6) singlet operator multiplied by a factor 

of m. This term is crucial since it cancels the negative potential energy for (j)^ introduced by 
O2. It is required by the central extension of the M = 2 supersymmetry algebra.'* 

Since the UV fixed point describes the conformal AA = 4 theory, whose coupling is a 
free parameter, one has the freedom to choose the UV coupling (or A) arbitrarily.^ We 
choose the 't Hooft coupling to be large, A S> 1, but finite as Nc — )• c«.^ The value of this 
coupling, together with the mass m, determines the intrinsic scale where the theory becomes 
strongly coupled, A ~ me"'*'^^/^ [22]. Notice, however, that for strong coupling the scale 
A is comparable to m and there is no large separation of scales. Consequently, one cannot 
decouple massive modes at the scale m while preserving the low energy spectrum at the scale 
A. 

The global symmetry group is SU{2)fi x U{1). The U{1) is an ordinary flavor symmetry 
under which $1 has charge +1, <I>2 has charge —1, and ^3 is neutral. The superpotential term 
mtr ($i$2) does not preserve the classical U{1)r symmetry that exists in the massless J\f = 4: 
theory. However, one can regard m as the expectation value of a background field and assign 

"^To see this, consider the identity {Qa, Qf} = ^a/S . To obtain a massive multiplet which has 4 helicity 
states (2 integer and 2 odd half-integer), the Bogolmolny bound must be saturated: Z^'^ = —Z^^ = 2m. This 
value for the antisymmetric matrix Z^^ forces half of the 2J\f fermionic "oscillators" (obtained from suitable 
linear combinations of the supercharges) to obey trivial anticommutation relations. This means that they 
cannot be used to create or annihilate additional helicity states. This, in turn, shortens a massive irreducible 
representation of the A/'-extended super- Poincare algebra from dimension 2^''^ to 2^, which is precisely what 
is needed to have a massive hypermultiplet. The action of and applied successively to O2 must yield 
the bosonic operator O3, and the anticommutation relation implies that an explicit factor of m must appear 
in O3. 

The dual operator analogous to O3 in the case of fundamental representation hypermultiplets is given 
explicitly in Appendix A of ref. [31]. Specifically, in eq. (A.l) of this reference, the operator Cm has an 
identical fermion bilinear, similar scalar trilinears, and similar mass terms for the hypermultiplet scalars. 
However, unlike O3, commutators do not appear in the scalar trilinears of Om- The reason is that their 
superpotential (in our notation) is tr ("l>i$3"l>2) rather than tr ("l>3[cE>i, c[>2]). We thank Andreas Karch for 
pointing this out to us. 

^One may also choose an arbitrary 9 angle, but this will not be important in our discussion. 

®This is necessary since the supergravity approximation to the string partition function is only valid for 
spacetime curvatures which are small when measured in string units. In our case, the ratio £a/L of the 
fundamental string length to the characteristic curvature radius L equals A~^'^*. 
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it an i?-charge of 2. This implies that the low energy effective action must be holomorphic in 
m. The vacuum manifold is an (A'^c— l)-dimensional complex space specified by adjoint scalar 
expectation values having the form (<I>i) = {^2) = and {^3) = diag(ai, . . . ^a^J (up to a 
gauge transformation), subject to the tracelessness constraint, Yli^i — ^- At a generic point 
in moduli space, Af = 2 supersymmetry is unbroken, the gauge group is Higgsed to U{1)^''~^, 
and the phase is Coulombic. A Wilsonian effective action for the low-energy Abelian the- 
ory may be generated by integrating out massive W bosons, their superpartners, and the 
adjoint hypermultiplet. Two derivative terms in the effective action are determined from 
a prepotential which depends holomorphically on the complexified gauge coupling constant 
r = 6/{2-k) + Ani/g'^, the mass m, and the eigenvalues {oj} [32]. The prepotential receives 
classical and quantum contributions. The perturbative correction is one-loop exact and may 
be determined by a matching calculation. Nonperturbative corrections become important 
in regions of moduli space where BPS states become light. Both monopole and W boson 
masses are proportional to differences of scalar eigenvalues, \ai—aj\, and vanish when eigen- 
values coincide. Away from points of eigenvalue degeneracy, nonperturbative corrections are 
suppressed and the low energy dynamics is well-described by the perturbative prepotential. 

In gauge/gravity duality, certain gravitational backgrounds may be interpreted as dual 
descriptions of AA = 4 SYM perturbed by relevant operators. In particular, solutions of 5D 
maximally supersymmetric M = 8 gauged supergravity which asymptotically approach AdS^ 
are dual to states of = 4 SYM with gauge group SU (Nc). If fields in addition to the metric 
are non- vanishing, then one can obtain non-conformal boundary theories with non-trivial 
renormalization group flow. Because 5D gauged supergravity is a consistent truncation of 
lOD type IIB supergravity on AdS^ x S*^, it is sufficient to solve the dimensionally-reduced 
problem (as opposed to the full lOD problem) [33-35] . It is possible to unambiguously "uplift" 
any 5D solution to a full lOD solution although, in practice, finding the uplift is nontrivial. 
The scalar and gravity part of the 5D action has the form 



Specific terms in £matter will be discussed below. Massless bosonic modes in lOD may be 
expanded in a set of Kaluza-Klein modes on the S^, with masses of order 1/L, where L is 
the radius of the five-sphere. There are 42 scalar fields and a complicated potential V which 
depends on 40 of them [33-35]. For a thorough description of the subsector we will study 
here, see refs. [23, 36]. The scalars fall into various representations of the 50(6) gauge group: 
there is a 20' representation with mass-squared = — 4/L^, a 10© 10 representation with 
= — 3/-L^, and dilaton and axion singlets which have = and do not contribute any 
potential energy. The scaling dimension A of a CFT operator Oa dual to a supergravity 
scalar field with mass-squared on AdS^ is given by A(A — 4) = M'^LP'. The 50(6) 
symmetry maps to the global i?-symmetry in the = 4 SYM theory. The 20' is dual to 
dimension 2 symmetric traceless combinations of the six AA = 4 scalars (i.e., the real and 
imaginary parts of (/>,), the 10 © 10 is dual to dimension 3 symmetric bilinears of the four 




(2.4) 
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TV = 4 Weyl fermions plus supersymmetric completions, and the dilaton-axion pair is dual to 
the M = 4 Lagrange density plus theta term.^ The scalar fields have < (but above the 
Breitenlohner-Freedman stability bound), which corresponds to A < 4. Therefore, a solution 
of 5D M = 8 gauged supergravity with non- vanishing profiles for any of the 40 nontrivial 
scalars corresponds to a relevant deformation of the boundary CFT. 

In our case, the deformations are encoded in two real supergravity scalars a and X: in 
terms of which the 5D matter Lagrange density equals 

-^matter = 3{daf + (Oxf + V{a, x), (2.5) 

with 



V{a, x) = r [-ie~"" - 5^'° cosh 2x + j-^ sinh^ 2x\ • (2.6) 

The dimensionful gauged supergravity coupling is if' = (2/L)^, and the 5D Newton's constant 
is Gs = 7rLV(2iVc). The field a has = -A/L'^ and is dual to the TV = 2* operator 
O2 in eq. (2.2). The second scalar field x has = — 3/L^ and is dual to the N = 2* 
operator in eq. (2.3). The value of the mass m in the field theory is determined by the 
asymptotic boundary conditions of the fields a and x^ which are not independent but related 
by supersymmetry. Notice that the identification of x with Os implies that the map between 
the field and the dual operator depends on the choice of boundary conditions/couplings in the 
field theory Lagrangian. This modification of the usual prescription may be a consequence of 
using a Kaluza-Klein compactification of ten-dimensional supergravity, we consider preserving 
supersymmetry on both sides as an indication that the identification is correct. 

A solution to 5D = 8 gauged supergravity with non- vanishing scalars a and x is known, 
and its full lOD uplift has been constructed [23]. The bulk 5D geometry has four-dimensional 
Poincare invariance. An analysis using a slowly-moving probe D3-brane indicates that the 
background branes form a locus around the origin of the probe's moduli space (the space where 
the probe's potential vanishes) where the probe's kinetic energy vanishes. This 'enhangon' 
geometry is dual to a special vacuum of the J\f = 2* Coulomb branch where extra dyonic states 
become massless [32, 37]. We will see that this solution plays an important role in the low 
temperature regime. 



3 Equilibrium thermodynamics 

Buchel, Deakin, Kerner, and Liu numerically solved the 5D Af = 8 gauged supergravity equa- 
tions for a background with three-dimensional rotational invariance and a regular black brane 
horizon [26].*^ They also verified that the lOD uplift solves the IIB supergravity equations 
of motion. This geometry is dual to a thermal equilibrium state of strongly coupled J\f = 2* 

''The boundary value of the 5D dilaton-axion pair is the UV marginal coupling r. The physical running 
coupling, dual to the lOD dilaton and axion, is a nontrivial function of r, m, and the energy scale [23]. 
*The zero temperature solutions were previously constructed analytically by Pilch and Warner [23]. 
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mIT Tim 
Figure 1. Left: Pressure divided by ^N^T^, as a function of m/T. Riglit: Low temperature behavior 
plotted as a function of T/m. The dotted line shows the extrapolation of a quadratic fit to the ten 
lowest-temperature data points. 



SYM with gauge group SU{Nc) in the large Nc limit. There is a single dimensionless pa- 
rameter, m/T, which may be freely adjusted. Thermodynamic quantities, when divided by 
appropriate powers of T, must be dimensionless functions of m/T. 

We have extended the analysis of ref. [26] to a wider range of m/T . The holographic 
representation of the entropy and free energy is summarized in Appendix A. As m/T increases, 
the differential equations one must solve become increasingly stiff, making the numerical 
exploration of very large values of m/T challenging. 

The pressure (equal to minus the free energy density) may be expressed in the form 

p='^NlT^f{m/T). (3.1) 

The scaling function / is plotted in figure 1 . The free energy is obtained from the renormalized 
Euclidean supergravity action (see Appendix A). The prefactors in eq. (3.1) have been chosen 
to equal the AA = 4 value of the pressure, so that /(O) = 1. A high temperature expansion 
yields / = 1 — 27r~^ r(|)^ (m/T)^ + 0[(m/T)^] [24]. Our numerical computations extend 
out to m/T ~ 33.3, at which point / ~ 0.1037. Throughout the range of temperatures 
< m/T < 33, our numerically-determined black brane solution describes a thermodynamic 
state which is at least locally stable. Both the pressure and the specific heat (discussed below) 
are positive, and we presume that our solution is describing the true equilibrium state. As 
shown in the right-hand plot of figure 1, extrapolation of our numerical results clearly suggest 
that / vanishes linearly as T/m — )• 0, implying that the free energy density F/V = 0{T^/m) 
at low temperature. The interpretation of this scaling behavior is discussed below. 
The entropy density may be expressed as 

S/V='^N^T^a{m/T), (3.2) 

where the prefactor equals the entropy density of AA = 4 SYM [38], so that a{0) = 1. The 
entropy may be obtained from the horizon area (A. 8) or from a thermodynamic derivative 
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Figure 2. The entropy density, divided by ^ N^T^, as a function of m/T. 
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Figure 3. The specific heat to entropy ratio, Cy/S — d{\n S)/d{lnT), as a function of m/T. 

of the free energy, 5 = —dF/dT, which imphes that cr = f — \y f ■ The resulting scahng 
function a is plotted in figure 2. For m/T <C 1, a high temperature expansion yields a = 
1 — 7r~^r(|)^ [m/T)'^ + 0[(m/T)^] [24]. As the temperature drops, a decreases and appears 
to approach zero linearly in T/m, implying that the entropy density S/V = 0{T'^/m) at low 
temperature. 

The specific heat at constant volume is given by 

'''' 

A plot of the specific heat to entropy ratio, Cy/S, is shown in figure 3, from which it is 
evident that the specific heat is always positive. The speed of sound (squared) in the plasma 
is inversely related to the specific heat, c^^^^^^ = Cv/S. Fi gure 3 shows that c^^^^^ equals 
1/3 in the high temperature (or massless) limit, as required from conformal invariance, with 
corrections quadratic in m. Explicitly, c^ound ~ ^ ~ i'''""'^ -'^(f)^ {m/T)'^ + 0[{m/T)^]. At low 
temperatures, the speed of sound (squared) evidently approaches 1/4. 
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Figure 4. The scaling function for ttie equation of state, w = (e — 3p) /{N^T^), as a function of m/T. 

The equation of state, relating the energy density e and pressure p to the temperature, 
may be expressed as 

e-3p = N^T'^ioim/T), (3.4) 

where 

This function is plotted in figure 4. Since e — 3p is the trace of the energy-momentum tensor, 
it vanishes in the conformal limit. For small m/T the scaling function uj rises quadratically, 
io = 57r-2 r(|)'^(m/r)2 + 0[{m/Tf] [24]. As shown in figure 4, uj peaks at m/T ^ 4.83 and 
decreases monotonically thereafter. The position of the peak corresponds to the temperature 
below which the plasma is cool enough that massive components of the hypermultiplet have a 
low probability of being thermally excited. The high-mass/low-temperature tail of to must be 
predominately due to the non-renormalizable (i.e., higher-derivative) interactions of the light 
degrees of freedom. The fact that m/T w 33 lies well out in the tail of uj, far from the peak, 
is evidence that our numerical results are probing behavior deep into the low temperature 
regime. 

The low temperature regime is dominated by the special enhangon solution mentioned 
in Sec. 2. This is shown by the plot in figure 5, which displays the evolution of scalars fields 
(x,e^") along the radial direction for solutions at different values of m/T, superimposed 
on the T = enhangon solution, whose explicit form can be found in Appendix A. As the 
temperature drops to zero, our solutions clearly approach the enhangon solution. As discussed 
in ref . [39] , the dual to the enhangon geometry may be interpreted as a flow from the A/" = 4 
SYM theory in the UV to a five-dimensional CFT in the IR. At low temperatures the effective 
five-dimensional theory dominates the dynamics, so this explains the behavior of the free 
energy, F/V ~ T^/m, the entropy, S/V ~ T'^/m, and the speed of sound, cl^^^^^ — >• 1/4, that 
we observe.^ 

^This is also consistent with the resuh of refs. [27, 28] that in the zero temperature hmit the bulk viscosity 
saturates the bound (:^/7)>2(| — c^), in agreement with ref. [40]. 
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Figure 5. Evolution of the scalar fields (x^e^") as the radial coordinate varies. All curves start at 
(0, 1), corresponding to the boundary where the scalars a and % vanish. Curves end at the point 
(x/uS^"'') corresponding to the value of the scalars at the horizon. The T = enhangon solution is 
shown as a dashed curve extending from (0, 1) to (oo,0). Finite temperature solutions are shown as 
solid lines for m/T « 2.23 (blue), 4.57 (red), 10.0 (green), 17.1 (purple), and 33.3 (orange). Our finite 
temperature solutions clearly approach the enhangon solution as m/T — >■ oo. 



This low-temperature behavior is presumably specific to the Nc ^ oo limit. Turning on an 
infinitesimal, non-zero temperature will deform the moduli space of degenerate vacuum states 
of = 2* theory into some non-trivial free-energy surface. One expects the free energy to be 
minimized at those points in moduli space with the maximal number of massless degrees of 
freedom. In = 4 SYM, a non-zero temperature destabilizes all points on moduli space other 
than the origin, where all scalar expectation values vanish and the entire SU{Nc) gauge group 
is "unbroken." Conformal invariance guarantees that the non-Abelian plasma phase, with 
0{N^) free energy, extends all the way down to T = 0. For the non-conformal J\f = 2 SYM 
theory, the low temperature behavior is more subtle. There is no point on the N = 2 moduli 
space where a non-Abelian gauge field is present in the low energy dynamics. Low energy 
degrees of freedom consist of ?7(1)^'^^^ gauge fields plus, at special points in moduli space, 
massless BPS monopoles or dyons [41-43]. One would expect low temperature equilibrium 
states to correspond to points where a maximal number of BPS dyons simultaneously become 
massless. This suggests that for sufficiently low temperatures there will be multiple degenerate 
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equilibrium states, related by SL{2, Z) electro-magnetic duality [42]. Hence, for finite Nc, we 
do expect a distinct low temperature phase to exist for some non-zero range of temperature. 
However, it is quite plausible that the transition temperature scales as A/N" for some a > 0. 
This would be consistent with our numerical results shown in figure 1, and is also consistent 
with the expected domain of validity of the low-energy description near the degenerating 
points studied in ref. [44]. Further study of the low temperature behavior of this theory at 
large but finite values of would be interesting. 

4 Screening masses 

Using holography, one may compute the thermal screening mass associated with a gauge 
invariant operator O by solving for the lowest eigenvalue of the linear operator which describes 
fluctuations in the corresponding dual field [46-49] . Let us remind the reader of the underlying 
logic. The long-distance behavior of the Euclidean two-point function of O determines the 
screening mass. If the Fourier transform of the correlator, viewed as a function of the spatial 
wavevector, is analytic in a strip of width k above the real axis, then the coordinate space 
correlator will fall with spatial separation |x— x'| at least as fast as g-'^l'^-'^'l. Hence, the 
screening mass k equals the distance from the real axis to the nearest singularity. If the 
Fourier transformed correlator is analytic except for poles on the imaginary axis, as will be 
the case in the large Nc theory under consideration, then the associated screening mass equals 
the magnitude of the pole closest to the origin. 

There is a well-defined prescription for calculating Green's functions using gauge/gravity 
duality. One must find solutions to the linearized equations of motion for the fiuctuations 
of the supergravity field dual to O. We will use combinations of fluctuation fields which 
are invariant under linearized diffeomorphisms. After a 4D Fourier transform, fluctuation 
fields depend on the radial coordinate r, the spatial wavevector k, and a discrete Matsubara 
frequency Un = 27rnT. We will limit attention to zero frequency fluctuations, since non- 
zero frequency fluctuations have shorter correlation lengths than the corresponding static 
fluctuations. 

The resulting linearized fluctuation equations of motion are second order ordinary dif- 
ferential equations (ODEs) in r. For some channels there is a single ODE involving just one 
fluctuation field; for other channels there is a coupled set of ODEs involving multiple fields. 
The ODE coefficients depend on the radial profiles of the scalar fields a and x plus the 5D 
metric components describing the background geometry, as well as the fiuctuation wavevector 
k which appears in the dimensionless combination (kL)^. 

For the sake of discussion, consider a gauge invariant supergravity field Z{r) that obeys a 
second order ODE. A general solution may be expressed as a linear combination of two inde- 
pendent solutions having simple behavior near the horizon, Z{r) = 0^°^ Z^°'^{r) +0^^^ Z||°''(r), 

^°For Af = 2 gauge theory, a detailed analysis of the temperature- induced deformation of moduli space has 
only been done for the case of an SU{2) gauge group [45], where it was argued that the discrete Z2 -R-symmetry 
is spontaneously broken at low temperature. 
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with coefficients c^°'' which will depend on the wavevector k. The horizon is a regular singular 
point of the fluctuation ODE and one solution, which we will choose to call Z^^, will have 
a logarithmic singularity at the horizon. The physical solutions in which we are interested 
should be regular at the horizon, and hence we must set cf™ = 0. The general solution may 
also be expressed as a linear combination of two solutions having simple behavior near the 
boundary (r oo), Z{r) = o^'^-^ Z^'^^{r) + o^^^ Z^^^{r). We may choose independent solu- 
tions such that Z^'^^ vanishes as r — t- oo while Zj^^ — )• 1 as r — >■ cx). Hence, the coefficient 
c^i^ equals the asymptotic value of the general solution Z. 

The asymptotic value o^^^ functions as the source for the (Fourier-transformed) gauge 
theory operator O, in accordance with the standard AdS/CFT interpretation. Hence, the 
Euclidean two-point function for O is the second functional derivative of the regulated on- 
shell supergravity action with respect to the value of Z at the r = oo boundary. The part 
of the on-shell action quadratic in fluctuations reduces to a boundary term of the form 
limr_>oo / d^fc J^(k)Z'(r)Z(r) plus contact terms which do not contain Z'{r), where the func- 
tion -F(k) depends on details of the action [50]. This leads to the result (OO) ~ c^'^^ /c^^^ , up 
to terms analytic in k'^ which can be removed by local counterterms.^^ Poles of the Green's 
function correspond to zeros of c^^^ . Consequently, to determine the screening mass one 
must find the least negative value of (kL)^ for which the linearized equation of motion for the 
fluctuation Z has a solution satisfying a Dirichlet condition at the boundary together with 
regularity at the horizon, cff"^ = c^^-^ = 0. 

We apply this procedure to evaluate the lightest screening masses for the energy- momentum 
tensor T^i, oi M = 2* SYM. This operator is dual to fluctuations of the 5D metric in A^ = 8 
gauged supergravity. We also obtain screening masses for the AA = 4 Lagrange density, 
vCa/'=4 = i tr F^^F^ji_y + • • • , dual to the 5D dilaton, as well as the Pontryagin density, 
iiE ■ B = \ tr F'^'F^u, which is dual to the 5D axion. 

In each Euclidean correlator the operators are spatially separated along the direction, 
so held theory coordinates transverse to this longitudinal direction are x^, x^, and r. Distinct 
symmetry channels are classified according to irreducible representations of the 0(2) group 
of rotations in the 1-2 plane, together with discrete eigenvalues for Euclidean time reversal 
(taking r — )■ — r), which we will denote as TZr, and charge conjugation C. (The eigenvalues 
of TZr and C will be denoted by Rr and C, respectively.) Irreducible representations of 0(2) 
are labeled by a non-negative integer angular momentum (or helicity) S plus, when 3 = 0, 
an eigenvalue for reflections in the 1-2 plane. We will use Ry to denote the eigenvalue of the 
reflection TZy taking — )• — x^. The resulting 3^^^ assignments for various operators and 
supergravity fields are conveniently tabulated in ref. [17]. 

Using coordinates x^ = (r, x^, x^, x^, z) with L = 1, the Euclidean black brane metric 
has the form 

ds^ = e2^(^) [B{z f dr^ + dx^] + . (4.1) 



^^The ratio c^'^^ /c^^^ is fixed by the required regularity of the solution at the horizon; this ratio does not 
depend on the actual value of Cji^. 
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We have switched to an inverted radial coordinate z = r^/r which places the horizon at z = 1 
and the boundary at z = 0. At the horizon, B[z) vanishes while A{z) remains finite and 
non-zero. Asymptotically, A{z) ~ — Inz and B[z) ~ 1 as z — )• 0. Appendix A discusses the 
background geometry in more detail. 

We add to this background a perturbation with arbitrary radial dependence multiplying a 
plane wave in the direction with imaginary wavenumber k. In other words, g^y = g^^ + h^^, 
where g^^ is the background metric (4.1) and 



^3 



V = V(^)e""" • (4.2) 

We similarly expand the two scalar fields present in the supergravity dual, writing a = Od + a 
and X = Xc\ + X? where the classical parts solve their respective Klein-Gordon equations [with 
potential (2.6)] and the perturbations, indicated by tildes, have the form 

a^5iz)e-'^-\ x^me-^""" ■ (4-3) 

Under 0(2) rotations, the various fluctuation components transform either as scalars (3 = 0), 
vectors {8 = 1), or rank-2 tensors {3 = 2), as shown in table 1. The 0(2) symmetry guarantees 
that the equations for fluctuations with differing helicities decouple. 

An infinitesimal diffeomorphism, 6x^ = S,^, produces a metric perturbation 6hf^i, = 
—^fiCu — ^uC^i, along with variations of the scalar fields given by 6a = —S,^d^aci and 
6x = —i^dnXcX- The covariant derivative V = d^S^y — T^^^a is taken with respect to 
the background metric. Of particular relevance will be residual diffeomorphisms which pre- 
serve the functional form of the metric perturbation, namely = £,^{z) e~'^^ . It is helpful to 
note that there exist linear combinations of metric perturbations (without radial derivatives) 
which are invariant under these residual diffeomorphisms. For example, hi2 and /loa are resid- 
ual diffeomorphism invariant perturbations in the tensor and vector channels, respectively. 
In the scalar channel, — /ioo + X^^=i K^ooZ-'^li) ^+ invariant under residual diffeomorphisms. 
We will fix this residual diffeomorphism invariance by imposing the "axial" gauge condition, 

/iM = 0, (4.4) 

at all points in space. 

The Einstein equations resulting from the action (2.4) are 

R^y = 4 {3df,a dya + d^x dvX) + |^(a, x) dtiv , (4-5) 
with the potential V{a,x) given by eq. (2.6). The scalar field equations are 

^ df,{V^g^''dya) = l:Va, -^d,{^gg^^''dyx) = \V^, (4.6) 



where here and henceforth we use subscripts to denote partial derivatives of the potential, 
Va = 5V/5a, = SV/5x, etc. Linearizing the Einstein equations about the background 
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3 


mode 


fields 


2 


tensor 


hab — ^h+5ab 


1 


vector 


hoa, ha3, haA 





scalar 


hoo, ^03, ^04, ^33) ^34, ^44, X 



Table 1. Classifications of supergravity fields under 0{2) spatial rotations. The index a = 1,2 
labels transverse spatial directions, 3 is the longitudinal spatial direction, and 4 the radial direction; 
h+ = hii + h22- 

metric g^^ produces the small fluctuation equation 



livh 



+ Mdf.Xci duX + d^x duXci) ■ 
The linearized scalar fluctuation equations resulting from eq. (4.6) are 



5^ iVa^i {d^X - h^'du Xci)] + \ d^h d^xci = I {V^x X + S) . 



5cl 
1 



(4.7) 

(4.8a) 
(4.8b) 



In these fluctuation equations, the background metric 5^'^ is used to raise and lower indices, 
compute covariant derivatives, and define the trace h = g^^ih^^. Scalar potentials are evalu- 
ated on the classical solutions ad and Xci- 

Not counting the dilaton and axion, there are two equations for the scalars a and x 
and fifteen equations for the components of the metric (which is a symmetric tensor). In 
the axial gauge (4.4), five of these equations become constraints and there are actually only 
five independent components of the metric, which split according to helicity into two for the 
tensor channel, two in the vector channel, and one in the scalar channel — which couples to 
the scalar fields a and x- 



Tensor channel 

The tensor mode involves the symmetric traceless perturbation hab — ^h^Sab- Its two inde- 
pendent components are h- = hn — /122 and /112, which obey identical second order ODEs. 
The resulting equation may be written more compactly if we let b{z) = lni?(z),^^ and use 
redefined fields H- = e~'^^h- and H12 = e~^'^/ii2. The result is 



H'L + {4.A' + h' + z^^)H'_ + K 



2 -2 -2A 

z e 



0, 



(4.9) 



^^This redefinition is not helpful for numerical calculations, since h{z) diverges at the horizon. Nevertheless, 
we use it here to make the equations more concise. 
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(and likewise for Hi2)- This is the same equation obeyed by the dilaton and axion — they 
are 'inert' scalars which have constant background values and obey source-free Klein-Gordon 
equations [51]. Linearly-independent near- horizon solutions have the form 

Hl{z)^l + ---, H^}{z)^{l + ---) + Hl{z)ln{l-z), (4.10) 

where ellipses denote ascending power series in z—1. As anticipated, the second solution has a 
logarithmic singularity at the horizon, and must be discarded. Since e~'^^ j z^ is regular at the 
horizon, it follows that the appropriate boundary conditions at the horizon are -ff-(l) = cq 
and H'_{1) = 0, with cq an arbitrary constant. Linearly-independent near-boundary solutions 
have the form 

HL{z) ^ z^{l + ■■■), Fi^(z)~ (! + •••) + (const.)i/L(z)lnz, (4.11) 

where these ellipses denote power series in z. The coefficient of the non-normalizable solution 
must be set to zero, so that H- vanishes as z^ for small z. The numerical analysis 
of the resulting boundary value problem is discussed in Appendix C. The lowest screening 
mass obtained from this analysis is plotted in figure 6 as a function of the ratio of the 
mass deformation to the temperature, m/T. As m/T — )• 0, we recover the AA = 4 value 
K w 3.4041 vrT. Extrapolating the low temperature data shown on the right in figure 6, we 
find k/{-kT) 3.248(5) as T/m O.^^ 

Vector channel 

In axial gauge the only independent vector fluctuation is /loa-^^ The Oa equations are fully 
decoupled. After defining Hoa = e^^^hoa, these equations take the form 

H'^, + {4A' -b' + z-^)H'oa + K^z-^e-^^Hoa = . (4.12) 

Near the horizon, linearly independent solutions are 

Hliz) ~ (z-l)2 (! + •••), Hliiz) ~ (1 + . . . ) + (const.) hI{z) In(l-z). (4.13) 

Both solutions are finite and once differentiable at z = 1. Consequently, one cannot select a 
regular solution by specifying a boundary value and first derivative precisely at the horizon. 

In performing the zero temperature extrapolation, we assume that the screening mass admits a Taylor 
expansion of the form '^^^QPn{T/m)"' . Truncating at quadratic order, we used a least-squares fit to the 
data starting from the smallest computed value of T/m. We looked for plateaus in plots of the coefficients 
Pn versus the number of fitted points, and found a very stable fit using 125 points corresponding to 0.0300 < 
T/m < 0.0315. The leading coefficient po ~ 3.2484. Deviations were scattered randomly around zero with 
displacements on the order of ±10~*. Our estimate of the uncertainty in the last digit is based on linear and 
cubic fits, which fit the data more poorly than a quadratic fit and predict values for po which under- and 
over-shoot 3.2484. 

^*The ha4 component is zero from our gauge choice and the (a4) equations become a constraint that can 
be solved as ha3{z) = Cate^^ for arbitrary constants Cai- Using the background equations of motion (A. 4), 
it is straightforward to show that this solution satisfies the second order (a3) equations. The undetermined 
constant in has refiects the fact that the axial gauge condition does not eliminate all diffeomorphism freedom. 
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Figure 6. Screening mass (in units of ttT) in the tensor channel, 2'^^^ = 2++, as a function of 
m/T (left) and, for low temperatures, replotted as a function of T/m (right). The dotted line is an 
extrapolation based on a quadratic fit to the data (see footnote 13). Fluctuations of the axion field, 
with ^'^-f^T — 0^~, and the dilaton, with g^^^^ = 0^^, obey exactly the same linear equation as the 
3 = 2 tensor mode, and consequently exhibit the same screening mass. 
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Figure 7. Screening mass (in units of ttT) for the vector channel, 3^^^ — 1"' , plotted as a function 
of m/T (left) and, for low temperatures, replotted as a function of T/m (right). The dotted line is an 
extrapolation based on a quadratic fit to the data. 



However, one may use the regular near-horizon asymptotic solution Hq^ to generate boundary 
conditions at some value of z close to, but slightly less than 1.^^ The near-boundary limit 
of the mode Hoa is identical to that of since b' ^ + 0{z^). Once again, we demand 
that HQa vanish like as z — >• 0"*" . Numerical implementation of this boundary value prob- 
lem is identical to that for the tensor mode. Our results for the resulting vector channel 
screening mass are plotted in figure 7. The limiting value for m/T = is k ~ 4.3215 vrT.^^ 
Extrapolating to zero temperature, we estimate that k/{ttT) — )• 3.999(9) as T/m — )• 0. 

^^Alternatively, one could remove the ambiguity by redefining Hoa{z) = {z—l)Hoa{z). This shifts the roots 
of the indicial equation down by one unit, from {2, 0} to {1, —1}, so that the regular solution now behaves as 
Hoa{z) ~ z — 1, giving boundary conditions Hoa{l) = and HQa{l) = 1, while the irregular solution diverges. 

^^This value is slightly smaller than the Af — i value quoted in ref. [17] which was obtained in ref. [49]. The 
discrepancy of ~ 0.0002 nT presumably reflects less accurate numerical work in these references. 
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Scalar channel 

In axial gauge, the scalar channel couples the metric perturbations /iqO) ^+ = ^ii + ^22) and 
/133 with the scalar perturbations 5 and x-^^ The mixing of these fluctuations at nonzero 
m/T is inevitable — the conformal Ward identity given in eq. (A. 31) involves the thermal 
averages of the operators O2 (dual to 5), O3 (dual to x)i and the trace T'^ (dual to a linear 
combination of /iqo, and /133). It is convenient to make the following field redefinitions. 
First, let /i+ = e^^H+ and /iqo = e^^^+^^i^oo- Second, define 

Z^H^^-\{l + b'/A')H+, a^a-\{a'JA!)H+, c ^ x - \{xJA') H+ . (4.14) 

These fields are invariant under infinitesimal diffeomorphisms preserving axial gauge. In 
Appendix C we derive the coupled system of ODEs for these gauge-invariant variables. The 
result is 

Z" +pZ' + qZ = faa + fcc, (4.15a) 
a"+paa' + qaa = gz{Z' + b'Z)+gcC, (4.15b) 
c" +PcC + qcC = hz{Z' + h' Z) + haU, (4.15c) 



where 



AA! + b' + z-^ -Ih'V/Q, fa = lh'[{?,A'+h')Vc, + 2Aa'^^V]/Q, (4.16a) 



p 


= AA! + b' + z- 


q 


= ^2^-2e-2A _ 


Pa 


= 4yl' + b' + Z-' 


Pc 


= AA! + b' + Z-' 


qa 


= K^z-^e-^^ - 


qc 


= K^z-^-^^ - 


ac 


— 6 ^<^X ^ 3 


ha 


- 1 r-2y I 4 



8 A'2 
3 



h'^V/Q , /, = I 6' [(3A'+6') + 8 Xd^^] /Q, (4.16b) 

gz = {\^V^ + ^^a'ciV) IQ, (4.16c) 
/iz = (^^% + ix'ci^) /Q, (4.16d) 



\z-^Vc,a - I [(6A'+6') a^iKv + 24 a'J y] /Q, (4.16e) 
h'^V^x - I [(6A'+6') x'ci^x + 8x'J^] /Q, (4.16f) 

[(3^'+6') a',,V^ + A'x'ciK. + 8 a^iXd^] /Q, (4.16g) 
[(3A'+6')XdK. + ^^a,,V^ + 24 a^ix^i^/] /Q , (4.16h) 

with Q = z2A'(3A' + 6').^^ 

It is instructive to first consider the massless limit, m = 0, corresponding to the AA = 4 
theory where we know that Od = Xci = 0. This implies that V = —\g^ = —3, Va = V^ = 0, 
Vaa = —6^^ = —24, y^-^ = —§5^ = —6, and ~ ^- source terms on the right 

hand sides of eqs. (4.15) vanish, leaving fully decoupled ODEs. The scalar fields a and c obey 



^'^The perturbation /103 may be found by integrating the first order (04) equation, leading to /io3(z) = 
Co3-B^e^^ for an arbitrary constant C03. Using the background equations of motion (A. 4), one may check that 
this solution also satisfies the second order (03) equation. 

^**Equations (4.f5)-(4.f6) are symmetric under the interchange ^/Za -f-)- x- The factor of reflects the 
normalization of the scalar kinetic terms in the action (2.5). 
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an equation of the form (/>" + [AA' + b' + 1/^)0' + (k^c-^^/z^ - M'^/z'^)(j) = 0, with = -4 
and —3, respectively. The leading behavior of linearly independent solutions is 

near horizon: a^{z) ~ IH , a"(z) ~ (H ) + a^{z) In(l-z), (4.17a) 

c\z)^ ! + ■■■, c^\z)^{l + ---) + c\z)\n{l-z), (4.17b) 

near boundary: a^{z) ~ z"^ 

(1 + •••), a^\z) ^ z'^{l + ---) + a\z)lnz, (4.17c) 
c^(z) ~ 2^(1 + •••), c"(z) ~z(l + •••) + (const. )c^(z) In z. (4.17d) 

The term —M'^jz^ in the ODE is unimportant near z = 1, so both scalars have the same near- 
horizon behavior. The metric perturbation field Z obeys Z" + (4A' + 6' + z"-*^ + 2+{zA'y^ )^' ~^ 
(k^z^^c"^"^ + 2^f'zA')^ ) ^ = 0- The linearly independent solutions are 



near horizon: Z^(z) ~ Z"(z) ~ ^ + Z^(z) In(l-z), (4.18a) 



near boundary: Z^(z) ~ z^(l + • • • ), ^"(z) ~ (1 + •••) + (const.) Z^(z) In z . (4.18b) 



For numerical evaluations (at arbitrary m/T)^ it is convenient to make field redefinitions 
such that the unwanted solution (or its first derivative) diverges at the horizon and boundary, 
while the physical solution remains regular. Examining eqs. (4.17) and (4.18), one sees that 
this can be accomplished by letting 

Z = z^{z-l)-'^Z , a = z'^a, c=z^c. (4.19) 

It is also convenient to define the vector 

W(z) = {Z', a', c', Z, a, c), (4.20) 

in terms of which eqs. (4.15) become a system of six first order ODEs. One can write this 
homogeneous system as a matrix equation of the form 

W' = L(z;k2)W, (4.21) 

where L is a 6 x 6 matrix involving the coefficients (4.16) and field redefinitions (4.19). For 
any choice of the vector W at the boundary, the solution W(z) at an arbitrary point in the 
bulk is linearly related, W(z) = ^(z; k^) W(0), where the transfer matrix $ depends on k? 
and equals the identity at the boundary, $(0; k^) = 1. Expanding the solutions in power 
series near the horizon and boundary, one can show that 

W(l) = (-2Zo, -2So, -2co, Zo, Sq, 2o), (4.22a) 
W(0) = {Zo,o, 0, 2o,o, 0, ao,o, 0). (4.22b) 

There are three undetermined coefficients at the horizon (Zq, oq, cq) and three at the bound- 
ary {Zofi, aofi, co,o)- In other words, at either endpoint there is a three-dimensional subspace 
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Figure 8. Screening masses (in units of ttT) for tlie scalar channel, 3f( ^ = 0^^, plotted as a function 
of m/T (left). For low temperatures, the lightest screening mass is replotted (right) as a function of 
T/m. The dotted line is an extrapolation based on a quadratic fit to the data. 



of solutions satisfying physical boundary conditions at that endpoint. Our task is to find 
those values of for which the boundary-to-horizon transfer matrix maps some vector in the 
physical subspace at the boundary into a vector lying in the physical subspace at the horizon. 
This may be accomplished by testing linear dependence. We choose basis vectors that span 
the physical subspaces, 

horizon: {(-2,0,0,1,0,0), (0,-2,0,0,1,0), (0,0,-2,0,0,1)}, (4.23a) 

boundary: {(1,0,0,0,0,0), (0,0,0,0,1,0), (0,0,1,0,0,0)}. (4.23b) 

Each basis vector defines a set of initial conditions which may be used to integrate the 
differential equation from one endpoint, either z = or 1, to a matching point in the bulk, 
z = z^. This yields a set of six vectors {W(z=i,)i=i_...^6}) where the first three are obtained by 
integrating from the horizon and the last three are obtained by integrating from the boundary. 
For generic values of ^ these six vectors will be linearly independent. But if, for some value 
of K^, a solution exists which satisfies the boundary conditions at both boundary and horizon, 
then the solution may be expressed either as a linear combination of {W(z*)j=i^2,3}, or as 
a linear combination of {W(z*)i=4,5,6}- In other words, for this value of k^, the six vectors 
{W(z=K)j=i^...^6} will not be linearly independent. If we regard these vectors as columns of a 
6x6 matrix, 

D= (W(z,)i, ... , W(z,)6) , (4.24) 

then det D must vanish for values of which correspond to screening masses. Consequently, 
the task reduces to finding the roots of detD as a function of ^^ . Appendix C describes our 
numerical procedures in somewhat more detail. 

Our results for the lowest screening masses are shown in figure 8. The operators dual 
to the scalar mode include the energy density Tqo as well as O2 and Oj,. When m/T = 0, 
the energy-momentum tensor and the scalar operators decouple. At this point, the values 
shown in figure 8 correspond, in ascending order, to the lowest screening mass of the ©2, 
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Too, and O3 operators in the undeformed A/" = 4 theory. The value of the second screening 
mass, K ~ 2.3361 vrT, associated with Tqo agrees with previous results [17]. For any nonzero 
value of m/T, these operators mix. Diagonalizing the matrix of correlators will, at long 
distances, yield decaying exponentials corresponding to the different screening masses shown 
on the left in figure 8. These curves may be regarded as the beginning of a tower of possible 
screening masses. The lightest screening mass is shown as a solid line, and its low temperature 
behavior is replotted as a function of T/m on the right side of figure 8. Extrapolating to zero 
temperature, we estimate that Hi/{-KT) — >• 1.620(2) as T/m — >• 0. 

5 Discussion 

Based on our results of screening masses in J\f = 2* SYM plasma, we make the following 
observations: 

• The smallest screening mass occurs in the Oljl^ channel for any value oi m/T. Fluctu- 
ations in the plasma with these quantum numbers are sourced by the energy density 
Too, as well as the scalar operators O2 and O3. These operators couple to the time-time 
component of the graviton and to quanta of the scalar fields a and x- The associ- 
ated screening mass is given by the solid curve in figure 8 which rises monotonically 
(in m/T) from 1.3731 vrT at high temperatures to 1.620 vrT at low temperatures. This 
screening mass should represent the mass gap, Kgap, whose inverse characterizes the 
longest possible correlations in the plasma. 

• The smallest screening mass in an Rr-oAd channel occurs in the 0^^ channel for any 
value oi m/T. It is shown in figure 6 and is our best candidate for the Debye mass. The 
associated supergravity mode is the 5D axion and its dual SYM operator is tv E ■ B. 

• Unlike the tensor and vector channels, the lightest screening mass in the scalar channel 
increases as the temperature decreases. The mass gap in the weakly coupled theory 
exhibits a similar behavior in the range of temperatures m ^ T S> A. The mass gap 
in units of the temperature depends quadratically on the coupling constant mgap ~ 
g'^{T)T, and for temperatures m ^ T the running of g{T) is the same as for an 
asymptotically free theory. 

• The screening mass extracted from the long-distance limit of the Too correlator jumps 
discontinuously when m/T goes from zero to an infinitesimal positive number. This 
happens for the following reason. The Euclidean two point function can be expressed as 
(Too(x)Too(O)) = e"'^"'''' |cnP, where {k^} may be thought of as excitation energies 
for the eigenstates {\n)} of the Hamiltonian (or transfer matrix) defined on X 
where the circumference of the compact spatial direction is /3 = 1/T, and c„ = (n|Too|0). 

^^To settle the question of whether this really is the gap, one should calculate the lightest screening masses 
in all possible symmetry channels. We have not considered SYM operators coupling to the NSNS and RR 
two-forms, since they have no simple 5D counterpart. 
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In the case of A/" = 4 SYM, the lowest energy with a nonvanishing ampHtude to couple 
to TOO is kq ~ 2.3361 ttT [17]. However, when a mass deformation is turned on, the 
SO{6)ji symmetry of the massless theory is explicitly broken to an SU{2)ji subgroup. 
Amplitudes c„ which vanish at m = 0, due to the S0{6) symmetry, may become non- 
vanishing at m 7^ 0, when the i?-symmetry is reduced. From figure 8 we learn that one 
of the corresponding energies is smaller than kq. Thus, the behavior of the correlator in 
the m/T — )• and |x| — )• oo limits depends on the order of limits. Examining the long- 
distance behavior first, for non-zero m/T, yields a Tqq screening mass which approaches 
1.373l7rT as m/T 0.^° 

Let us comment on similarities and differences between = 4 SYM, N = 2* SYM, and 
QCD plasmas: 

• In QCD, the mass gap and the Debye mass are found in the 0]jl~'' and channels, 
respectively. This is also true for our candidate mass gap and Debye mass in N = 2* 
SYM. And it is true for AA = 4 SYM with an important caveat: one must limit the 
comparison of operators to the -R-singlet sector. It is known that the screening mass for 
an i?-current in A/" = 4 SYM is less than that for the energy density [18]. However, the 
i?-current transforms nontrivially under the SO{Q)ii of AA = 4 SYM. QCD, of course, 
does not have such a symmetry, which is why the authors of ref. [17] chose to restrict 
comparisons to the i2-singlet sector. 

• In AA = 4 SYM, conformal invariance implies that screening masses are always propor- 
tional to the temperature. In pure Yang-Mills theory at finite temperature, screening 
masses divided by the temperature are nearly constant for 1.5Tc < T < 4Tc [19, 20]. 
Therefore, maximally supersymmetric and non-supersymmetric non-Abelian plasmas 
are likely to be most similar in this temperature window. In N = 2* SYM, there are 
two regimes of m/T where screening masses scale linearly with temperature: very high 
temperature, m/T <^ 1, where the mass deformation is negligible, and asymptotically 
low temperature, m/T — )• oo. These plateaus are clearly visible in, for example, figure 8. 
One might expect the low temperature J\f = 2* plasma to be most similar to QCD plasma 
at T ~ 2Tc, since the heavy M = 2 matter decouples in the low temperature regime of 
M = 2* theory. 

In table 2, we show screening masses in various symmetry channels, divided by vrT, in 
QCD [SU(3), A^f = 2, T « 2T^] [21], M = 2* SYM in the large m/T limit,^! and AA = 4 
SYM. As one sees from these results, screening masses in QCD at T ~ 2Tc are not small 
compared to vrT. This is a clear sign that QCD plasma in this regime is not weakly 
coupled. Screening masses in A/" = 4 SYM are roughly twice as large as in A'^f = 2 QCD 
at T ~ 2Tc. Larger screening masses (or shorter screening lengths) in AA = 4, relative to 
QCD, is to be expected since there are many more fields contributing to screening in 

■^"We thank Andreas Karch for discussions about this point. 

^'^More precisely, we show data from our largest mass value, m/T ~ 33.3. 
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qCRr 


N{ = 2 QCD 


M = 2* SYM 


AA = 4 SYM 


or 


1.25(2) 


1.62 


2.34 


ol- 


1.80(4) 


3.25 


3.40 


1+- 


2.88(12) 


4.01 


4.32 


2++ 


2.56(7) 


3.25 


3.40 



Table 2. Screening masses in selected symmetry channels, in imits of ttT, in QCD (A^f = 2, T w 2Tc), 
the large mass regime of TV = 2* SYM {m/T « 33.3), and A/" = 4 SYM. 

= 4 theory. As seen in table 2, screening masses in the mass-deformed M = 2* theory, 
at large m/T, are smaller than in AA = 4, and closer to the QCD values. The decrease 
in screening masses, relative to AA = 4, is relatively modest except for the 0^ scalar 
channel, where the change is substantial. 

We remind the reader that the Ol^ and 2~^~^ screening masses are necessarily identical 
in both AA = 4 and N = 2* SYM. This reflects the fact that in both these theories the 
axion obeys the same equation of motion as the transverse components of the graviton. 
A holographic theory in which these two supergravity modes obey distinct equations 
will be needed to more closely model QCD. 

• It is also interesting to compare ratios of screening masses in different symmetry chan- 
nels, instead of their absolute values. Figure 9 plots screening mass ratios, relative to 
the mass gap (or screening mass in the 0^"*" channel) as a function of m/T in M = 2* 
SYM. Table 3 shows the ratios of screening masses relative to the mass gap for QCD 
[Nf = 2, T 2Tc] [21], J\f = 2* SYM in the large m/T limit, and AA = 4 SYM. The ratio 
of the 01*^" and Olj]^ screening masses is virtually identical in QCD and A/" = 4 SYM. 
This near-perfect agreement is surely fortuitous, and in this particular case the M = 2* 
ratio deviates much farther from the QCD value. But, as shown in table 3, in the 
other symmetry channels for which we have results, the J\f = 2* screening mass ratios 
are much closer to the QCD values than the J\f = 4: ratios. 

Overall, considering both the absolute values and the ratios of screening masses in 
various symmetry channels, it seems fair to regard J\f = 2* SYM as a better model of a 
QCD plasma than AA = 4 SYM. 

• In pure Yang-Mills theory, screening masses divided by the temperature drop precipi- 
tously (at least in certain S = channels) when the temperature falls below 1.5rc and 
approaches Tc [20] . Neither AA = 4 nor M = 2* SYM have thermal phase transitions at 
a non-zero temperature. (For the former this is certain, and for the latter we find no 
evidence.) Thus, it is not sensible to ask how well either supersymmetric theory models 
QCD near Tc. However, as in QCD, screening masses of M = 2* SYM change substan- 
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Figure 9. Ratios of A/" = 2* screening masses (in indicated channels) relative to the mass gap (or 
screening mass in the 0^"*" channel), as a function of m/T . 



qCRr 


N{ = 2 QCD 


M = 2* SYM 


AA = 4 SYM 


ot- 


1.45(2) 


2.01 


1.46 




2.31(10) 


2.48 


1.85 


2++ 


2.05(6) 


2.01 


1.46 



Table 3. Ratios of screening masses in the indicated symmetry channel to the mass gap (or screening 
mass in the Oj+ channel), in QCD {Ni = 2, T « 2Tc), the large mass regime of Af=2* SYM 
(m/T w 33.3), and7V = 4 SYM. 

tially over a relatively small temperature range. From figures 6-9, it is clear that the 
high temperature plateau is limited to m/T < 1 and the low temperature plateau is 
roughly m/T > 10 — significant variation is confined to the band m/T k, few. The 
point m/T ~ 4.83, where the trace anomaly deviates maximally from zero (c.f. figure 4), 
lies in the middle of this band. 
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A Properties of the gravity dual 

In this and the following appendices, we present details of our analytical and numerical 
calculations at a level which would allow someone, with some effort, to reproduce the results 
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presented above. Further details, including explicit expressions for the coefficients of series 
expansions and more explicit derivations, may be obtained online [30]. 

This appendix summarizes the 5D J\f = 8 supergravity solution dual to finite temperature 
large- A'^c M = 2* SYM. Most of the discussion is taken directly from the original works [24-26]. 
We collect the results in one place for two reasons: to explain, in detail, how our calculations 
were performed, and to benefit readers wishing to do further studies. We used Mathematica 
7.0 for both analytical and numerical computations. 

Equations of motion 

We take the following ansatz for the 5D black brane metric, 

ds^ = e^Mr) ^_B{rfdt^ + dx^) + dr^. (A.l) 

Without loss of generality we have chosen coordinates such that grr = 1 everywhere. The 
metric is invariant under 3d spatial rotations, as well as time and spatial translations. The 
spacetime is defined for values of the radial coordinate r/j < r < oo. The lower end is an 
event horizon with the property that -B(r/j) = 0, and A{rh) = A^, a finite number. Infinity 
is a boundary where the metric approaches that of AdS^. This requires limr_s>oo B{r) = 1 
and A{r) to asymptotically approach r /L."^"^ The scalar fields are chosen to be functions of 
r only. They should be finite at the horizon: a{rh) = ah and x{'''h) = Xh^ and vanish at 
infinity: lim^^oo Oi{r) = \\mr^oo x{^) = 0. We can shift r and compensate the change in the 
asymptotic form of the metric by rescaling the spacetime coordinates. Let us then choose 
r/i = without loss of generality. The radial coordinate r is not well-suited to numerical 
calculations since the boundary is at infinity. Let us switch to a new, dimensionless radial 
coordinate z = e~'^l^ This places the horizon at 2 = 1 and the boundary at z = 0. In 
this new coordinate system, A{z) ~ — In 2: as z — )• 0. Since the logarithmic divergence is 
troublesome for numerical work it is helpful to extract this asymptotic behavior by defining 

= -lnz + I(z). (A.2) 

Note that this field redefinition does not significantly alter the behavior of the warp factor A 
near the horizon. In particular, A{X) = A^. The metric now reads 

p2A{z) J 2 

ds^ = ^-^{-B{zfdt^ + dx2) + l2 ^. (A.3) 

^^The asymptotic behavior of A{r) could also he r/L + A^, where is an arbitrary constant. However, 
Aoo can be set to zero by rescaling the time coordinate. 

■^^Some derivations are more transparent using the r coordinate, whereas technical calculations are better 
handled using the z coordinate. We will freely switch back and forth between r and z in the following. 
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The Einstein equations derived from the action (2.4) and metric (A. 3) are 



B' = Cz^e-^^, (A.4a) 

3z2 



[A' - \lzf + - \lz)B'lB - a'^ - + ^ = 0, (A.4b) 



a" +iAA! + B'/B-^/z)a' -'^^=Q, (A.4c) 

^' + (4l' + B'/B - 3/z)x' - 5 = 0, (A.4d) 

where a prime denotes d/dz and C is a constant that will be determined later. The system 

(A. 4) is sixth order. Note that factors of appearing in the system and in the gauged 

supergravity coupling present in the scalar potential cancel. Therefore, it is natural to set 
L = 1 and ^2 = 4.25 

Near-horizon solution: temperature and entropy 

To have a regular black brane horizon, near z = 1, we try a power series solution of the form 

B{z) = -Bh{z-l) [l + 5^ Sn(^-l)"] , (A.5a) 

n>l 

A{z) = Ah + Y,M^-'^T, (A.5b) 

n>l 

a{z) = ah + ^ an{z-lT, (A. 5c) 

n>l 

x{z) = Xh + Y.Xn{z-lT. (A.5d) 

ra>l 

Substituting these into system (A. 4) (first taking the logarithm of eq. (A. 4a) then differen- 
tiating with respect to z), we can solve the equations order- by-order in z — 1. We solve up 
to sixth order. The asymptotic values A^, B^, ah, and Xh are undetermined, and will be 
referred to as "horizon data." All higher coefficients {Bn, An, an, Xn} are fully determined by 
the horizon data. At the horizon, the fields and their slopes are finite, as desired, 

B{1) = 0, B'{1) = -Bh, a(l) = ah, a{l) = 0, (A.6a) 

A{l) = Ah, A'{1) = 1, x{l) = Xh, X(1) = 0. (A.6b) 

The Hawking temperature of the black brane may be extracted from the metric near the 
horizon using the usual analytical continuation to Euclidean signature with a compact time 



■^*By making the clever choice of radial coordinate t{z) = 1 — B{z), the authors of ref. [26] reduce the system 
to a set of three second order ODEs which does not explicitly contain B or its derivative. We decided not to 
use the t coordinate because it leads to near-horizon and near-boundary expansions with fractional powers of 
t; we find the series solutions simpler to understand in the z coordinate. 

^^This is different from the convention used in ref. [26] where L = 2 and = 1. 
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direction, and demanding that the solution is regular at the horizon, where the time circle 
collapses to zero size. This fixes the period of the temporal circle to be 

In the AdS/CFT correspondence, we identify the black brane temperature T = 1//3 with the 
temperature in the quantum field theory. 

The Bekenstein-Hawking entropy of the black brane is proportional to the volume of its 
three-dimensional horizon divided by Newton's constant in five dimensions, S = Vhorizon/ (4G5). 



The horizon volume is ^horizon = /r3 d^x det g^^^, where g^^^ = e^^'^Sab is the induced metric 
on the hypersurface z = 1 and fixed r. Therefore, Vhorizon = Ve^"^'*, where V is the (formally 
infinite) volume of x-space. The entropy density is 

Writing this as in eq. (3.2), we identify a = (2/3^)^. Positivity of entropy implies that > 0. 
We can fix the unknown constant C in eq. (A. 4a) using the near-horizon solution. Plugging 
eqs. (A. 5a) and (A. 5b) into eq. (A. 4a), then evaluating at z = 1, we obtain C = —e^^'^B^- 

Near-boundary solution: gauge-gravity dictionary 

Near 2 = 0, solutions may be expanded in an asymptotic series solution of the form 

00 n 

Biz) = l + ^5^i?2n+4,fc^'"+'(lnz)^ (A.9a) 

n=0 fe=0 
00 n+1 

^(^) = EE^2„+2,fc^'"+'(lnz)^ (A.9b) 

n=0 k=0 
00 n+1 

«(^) = EE"2n,fc^'"+'(ln^)^ (A.9c) 

n=0 fe=0 

00 n 

X(^) = EEx2n,fc^'"+ninz)^ (A.9d) 

n=0 k=0 

The leading (small z) behavior of the fields B and A are dictated by the requirement that, at 
the boundary, the metric approaches that of AdS^. Moreover, the two scalar fields must vanish 
at the boundary. The rate at which they vanish (i.e., the leading exponent) is determined 
by linearizing the scalar equations of motion at z = and solving the characteristic equation 
for an ansatz of the form z^.^^ The presence of Inz in the leading series solution comes from 



■^^The scalar potential is, to quadratic order in the small fields, —L^V « 3 + 12a^ + 3x^- Therefore, the 
linearized scalar equations are a " - (3/2)a' + (4/z^)a = and x" - (3/z)x' + (3/z^)x = 0. The first equation 
has p = 2 as a double root. Therefore, a pair of linearly independent solutions is and Inz. The second 
equation has p = {3, 1}, so a pair of linearly independent solutions is z'^ and z + cz^lnz, with c a nonzero 
coefficient determined by the equations of motion. 
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the fact that z = is a regular singular point and the roots of the characteristic equation 
differ by an integer. Higher powers of Inz arise from nonlinearities in the system (A. 4). The 
unknown coefficients are found by plugging the ansatz into the system (A. 4) and solving 
order-by-order in z°'{lnz)^. We carried this through up to order a = 8. Some coefficients 
vanish, ^2,1 = -Be,! = 0, but generically all terms in the series (A. 9) are non-zero. Observe 
that eq. (A. 4a) determines the i?„.fc in terms of the An^k, and eq. (A. 4b) determines the A^^k 
in terms of the an,k and Xn,k- The first order equations (A. 4a) and (A. 4b) should be thought 
of as constraints. The two scalar equations determine all higher an.k and Xn,k in terms of the 
first couple of coefficients. All higher coefficients {Bn^k, A^^k, ctn,k, Xn,k} are fully determined 
by the "boundary data": i?4,0) ao,i! oo,0) Xo,0! and X2,o- At the boundary, the fields and their 
slopes are 

B(0) = 1, B'(0) = 0, a(0) = 0, a'(0) = 0, 

^ (A.IO) 
A{0) = 0, ^'(0) = 0, x(0) = 0, x'(0) = Xo,o. 

We can fix the constant C in eq. (A. 4a) using the near-boundary solution. In fact, we may 
evaluate this equation at any value of z and the constant C must be the same — doing this 
at the horizon and boundary provides a powerful constraint. Plugging eqs. (A. 9a) and (A. 9b) 
into eq. (A. 4a), dividing by z^, then evaluating at z = 0, yields C = ^B^^. We now have a 
nontrivial equation that relates boundary and horizon data, 

454,0 = -e^'^'^Bh. (A.ll) 

To summarize, the near-boundary solutions are 

B{z) = l + BifiZ^ + 0{z^), (A.12a) 

A{z) = A2,oz^ + 0{z^ln^ z), (A.12b) 

a{z) = z^ [q!o,i In z + ao,o + 0{z^ In^ z)] , (A.12c) 

X{z) = z [xo,o + z^lxlo In z + X2,o) + 0{z^ In^ z)] . (A.12d) 

Generally, in the AdS/CFT correspondence, the leading coefficient in a near-boundary expan- 
sion corresponds to the source for the dual operator and the subleading coefficient corresponds 
to the expectation value of that operator. The operators dual to a and x are O2 and Os, 
respectively. Hence ao,i will be proportional to and qo,o is related to the thermal expec- 
tation {02)- Likewise, xo,o will be proportional to m and X2,o related to (O3). 

The precise correspondence relies on a matching calculation at zero temperature. At 
T = the background metric must be Poincare-invariant. This forces B = 1. The remaining 
fields A, a, and x obey a system of 1st order ODEs referred to as the supersymmetric flow 
equations [23], whose solution is 

= A:e2"/sinh2x, (A.13a) 
e^" = cosh2x + (7 + lntanhx)sinh2(2x), (A. 13b) 
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Xo.o 

Figure 10. m/T as a function of xo.o- 

where k and 7 are constants of integration. Plugging in the asymptotic formulas for the fields 
fixes the undetermined coefficients. At leading order, 

"0,1 = gXo.O' Xo,o = k/2. (A.14) 
At the first subleading order one finds 



ao,o 



5X0.0(1 + 27 + 2 In xo,o), X2,o = 2ao,oXo,o - ixno- (A.15) 



The constant 7 parameterizes a family of distinct solutions. Solutions with 7 < correspond 
to points on the Coulomb branch of = 2* SYM (where the scalars (pi and 02 are massive and 
03 gets an expectation value). Solutions with 7 > are unphysical. The expressions given 
in eq. (A.14) for the leading series coefficients should not change at finite temperature since 
they are bare couplings in a Lagrange density. As such, the constant k must equal mL times 
a pure number. Eq. (A.14) must also hold for the black brane solution. However, eq. (A.15) 
holds only at zero temperature. For the case 7 = 0, a probe D3-brane computation shows 
that the supergravity solution may be understood as a Coulomb branch vacuum in which the 
background D3-branes form a linear enhangon singularity. A computation of the size of the 
enhangon shows that k = mL. Thus, 

mL (mL)^ 
Xo,o = ao,i = — ^ — • (A.16) 

We choose to regard xo,o as an independent variable on which the other supergravity 
coefficients depend. In other words, the boundary and horizon data are all functions of Xo.o^ 

-B4,o(xo,o), ao,o(xo,o), X2,o(xo,o), ^/i(xo,o), a/i(xo,o), Xhixofl)- (A.17) 

There are six independent data — ao,i is fixed by eq. (A.16) and is fixed by eq. (A. 11). 
Solving for L in eq. (A. 7), the map from xo,o in sq. (A.16) to m/T is 

, 47rxo,o 1o^ 
mT= . A.18 
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A plot of this map is shown in figure 10. The mapping is nonhnear and monotonic, with larger 
values of xo,o corresponding to larger values of m/T. Since Ah{0) and 5^(0) are nonzero (see 
Appendix B), m/T = when xo,o = 0. 

Holographic renormalization, thermodynamics and expectation values 

The Euclidean supergravity action is regulated and renormalized following the approach of 
ref. [25]. Proper renormalization is needed to define thermal expectation values correctly. 
We start by cutting off the boundary at the hypersurface given by the equation r = tq. To 
the bulk Euclidean action S'buik given by (minus) eq. (2.4), we add the Gibbons-Hawking 
boundary term 

5bdy = -^^^/ d^xVhK, (A. 19) 



r=ro 



where h is the determinant of the induced metric h^i, on the hypersurface, and K is the trace 
of the extrinsic curvature. The sum of these actions will diverge in the limit that rg — )■ oo. 
To cure this one must add counterterms at the boundary. When the boundary is fiat the 
counterterm action has the general form 



'-'ctr — 



2 

2^^S( C1+C3 a+C4X^+C5 a^+C6 ax^+cg hcio x"^ Inx+cn x' 

\ Inx 



(A.20) 



r=ro 



where x = ^/goo(roj = e'^^^°-'-B(ro) parameterizes the location of the boundary. Note that 
the counterterm action density is completely local and that only even powers of x appear 
(this is necessitated by the fact that only even powers of z appear in the near-boundary series 
for A, B, and a). By tuning the coefficients Cj=i^...^ii appropriately, the regulated action 
Sreg = Shuik + 'S'bdy + 'S'ctr may remain finite upon taking the limit tq — 00. 
The on-shell bulk action, with the Gibbons-Hawking term evaluates to 



'-'bulk 



dr 



BV 

= £^A-S, (A.21) 



where, using eqs. (A. 7) and (A. 8), the term obtained from the lower limit is just the entropy 
5. The quantity A diverges in the limit ro — )• 00. Switching to the z radial coordinate (so 
that the hypersurface sits at z = zq = e~'^''^^), we see that the leading singular behavior of A 
arises from e^^^^''^ ~ z^'^. The counterterm coefficients are chosen to remove the sensitivity 
of A on the cutoff in the limit ~^ 0. A straightforward expansion in powers of Zq^ and 
In zo leads to 

Cl = ^, C3 = 0, C4 = C5 = f , = 0, = Cio = -3^. (A.22) 

The coefficient cn multiplies a finite counterterm and corresponds to a shift in the subtraction 
scheme. It may be fixed unambiguously as follows. The renormalized supergravity action is 
identified with the gauge theory free energy divided by the temperature, ^ren = lim^g^o •S'reg = 
/3F. Using F = E — TS, the energy is given hy E = lim2,)_j.o A. At zero temperature. 
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the energy must vanish since the ground state is supersymmetric. For A to vanish in this 
hmit, cii must equal Ultimately, 



lim A = -x^ifi - fxo,o[x2,o - (2ao,oXo,o - 5X0,0)]- 

ZQ 



(A.23) 

Converting from 1/L to T, and using eq. (A. 11), the free energy density can be written 



as 



F/V 



B. 



4,0 



14 2 

X2,o Xofi + 3X0,0 ~ 2qo,o Xo,o 



(A.24) 



The overall factor —^N^T'^ is the free energy density of A/" = 4 SYM. The remainder of the 
expression constitutes the scaling function / in eq. (3.1). 

In ideal fluid hydrodynamics, the energy- momentum tensor T^'^ can be expressed in terms 
of the local fluid four-velocity n^(x) as T^"^ = (e + p)u^u'^ + pg^^'^, where the energy density 
e = E/V and the pressure p = —dF/dV. The trace is given by 



-(n) = e-3p=(ii; + 3F)/V 



-64,0 ^ Bh 



2 \3 



) (x2,oXo,o + gXo,o - 2ao,oXo,o) • 



(A.25) 



In a conformal theory, e = 3p. The factor following N^T'^ is the scaling function u in eq. (3.4). 

We now describe how to obtain supergravity formulas for the thermal expectation values 
of the operators O2 and O3 using methods described in refs. [52, 53]. These expectation values 
are given by functional derivatives of the renormalized supergravity action with respect to 
the couplings and m (which, for this discussion, should be thought of as independent 
couplings) , 



{02) = my 



_l5Sr, 



lim 



1 6Sreg Sa^ro) 

1 6Srcg 6x{ro) 
5m ro-5>oo f3V 5x{'fo) 5m 



lim 



(A.26a) 
(A.26b) 



Let us evaluate the first functional derivative in each of the above expressions. Computing 
the scalar functional derivatives for Sctr is straightforward with eq. (A. 20), so let us focus on 
•Sbuik- We have 



47rG, 
1 

4^ 
1 



'S'bulk = — , J cfix ^/g(^\R + /^matter 

d^x ^ 5^2 + \ABIB + 2i + \BIB + Zo? + + V{a, x) 
d^x^\\{k + K^^'K^^) + I (2^2 + ABjB) + -id? + + V{a, x) 



47rG5 
1 

4^ 



d^x^\\{-K^ + K^^'K^y) + + 2x^ 



+ 



d^x^gK 



r=rQ 

(A.27) 
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The functional derivatives of the action are to be evaluated at tq. They should be viewed as 
the canonical momenta conjugate to q and x- One finds, 



— e B[6a + zcsa + 2cs 



r=ro 



(5a(ro) 47rG5 V Inx 
<5Sreg _ /3V g4A^/2x + 2c4X + 4ciox'lnx + 4cnx' 



dxiro) 47rG5 
Continuing with the calculation, write 



(A.28a) 
(A.28b) 



r=ro 



6a(zr)) 6ar)i6a(zQ) r2 nr, ^/ 9, 9 ^ 

= -S^^ = i.S[ln.„ + 0(.§W.„)], (A.29a) 

-1^ = ^^= 2^o[l + 0{z,lnzo)]. (A.29b) 

Putting it all together and sending the cutoff to zero (it is convenient to do this in the z 
coordinate) yields 

iV2 i\r2 
(^2) = ao,o, (O3) = (X2,0 + hlo)- (A.30) 

As expected, the subleading coefficients of the scalar fields in 5D supergravity are directly 
related to the expectation values of the 4d gauge theory operators to which the leading 
coefficients couple. Lastly, one may check that the expectation values we have computed 
satisfy the conformal Ward identity, 

-(T';) = 2m2(02>+m(03>. (A.31) 

At zero temperature, one may insert the value of X2,o found in eq. (A. 15) and confirm that 
hmT->o (K) = 0. 

Endpoints of thermal scalar flows 

Figure 5 showed the trend that the thermal solutions for the scalars a and x became succes- 
sively better approximations to the zero temperature enhangon solution as m/T — t- cxd. In 
figure 11 we provide further evidence for this assertion by plotting the residuals to eq. (A. 13), 
evaluated at the horizon, for our numerical backgrounds. On the left is a plot of e^*^ — 
2xo,o e^^^ I sinh2x/i, and on the right a plot of e^"'' — [cosh2x?i + (7 + Intanhx/i) sinh^(2x;h)] 
for 7 = 0. In each plot the deviation from zero vanishes as we approach larger values of Xo,0) 
which corresponds to lower temperatures. 

First law of thermodynamics 

The first law of thermodynamics, dE = TdS, is equivalent to the statement that dF/dT = —S. 
Let us express this constraint in terms of supergravity parameters. For temporary convenience 
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Figure 11. Deviation of endpoints of tlicrmal flows from supersymmetric flow equations. Left: a plot 
of e^*^ — 2x0,06^"''/ sinh2x/i vs. Xo.o- This measures how well eq. (A. 13a) is satisfied at the horizon. 
Right: a plot of e^"'* — [cosh2x/i + (7 + Intanhx/j) sinh^(2x/i)] vs. xo.o- This measures how well 
eq. (A. 13b) is satisfied at the horizon. 



define Vl = X2,o Xo,o + 5X0 ~ 2ao,o Xo o- Differentiating the formula for the free energy density 
in eq. (A. 24) gives 

The entire term inside the braces must equal 1. Note that the temperature derivative acts 
on objects which are functions of xo,o- Using the chain rule, T-^ = T 3^^, we need 

T ^ = [A'j^ + B'JB, - l/xo,o] ■ (A.33) 

This relation comes from differentiating eq. (A. 18), with primes denoting d/dxofl- After 
straightforward manipulation, the first law condition becomes^^ 

5X2,0X0,0 - ix2,o - "0,0X0,0 + 2ao,oXo,o + |^4,o - ^^^4^ = 0. (A.34) 

This is a 1st order ODE with respect to the parameter xo,0) involving four of the six inde- 
pendent horizon and boundary data. Therefore, it is a nontrivial constraint which serves as 
a robust check on our numerics [26]. 

The left hand side of eq. (A.34) is plotted in figure 12. While no numerical computation 
with finite precision numbers can ever give a result which is exactly zero, the scale on our 
plot indicates that the first law is satisfied accurately to within 10~^ for the largest xo,o (or 
lowest temperatures) explored. The deviation from zero is due to discretization effects.^® 

^'^This simple form is based on the assumption that the relation between T and xoo is monotonic. We have 
verified numerically that B'f^/Bh + A'f^ 7^ l/xo,o- Notice that dxofi/dT is always nonzero since Bh is always 
positive; this follows from the positivity of the entropy. 

^*To demonstrate that the deviation is a discretization effect, we tried using a better approximation for the 
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B Numerical procedure 

We used a shooting technique to determine the independent boundary and horizon data (A. 17) 
as functions of Xo,o- The algorithm is as follows. For a given value of Xo,0) start with a trial 
set of independent boundary and horizon data, 

X(xo,o) = (^4,0, "0,0, X2,o, Ah, ah, Xh){xofi)- (B.l) 

The remaining parameters are fixed by ao,i = |xoo ^^'^ = —454^06"^"^''- One may now 
construct a near-horizon series solution using eq. (A. 5). Use this series to evaluate 

\{z) = {B,A,a,a',x,x){z) (B.2) 

at some Zmax close to the horizon. The output is some vector V(2;jnax)|series which provides 
initial conditions for the equations of motion given in eq. (A. 4). Integrate this system from 
Zinax down to a point in the middle of the bulk. Evaluating the fields and their derivatives 
at this point produces V(2;=i,)|hor-s>buik- Now repeat this process from the other direction. Con- 
struct the near-boundary series solution using eq. (A. 9) and use it to evaluate V(zniin) I series, 
where z^[^ is close to the boundary. This provides the initial conditions needed to integrate 
the system from Zmin up to z*. Evaluating the solution at the last point gives V(z*)|bdy-s>buik- 
Finally, take the difference 

M - V(z.)|,,^^,^j, - V(z,)|,^^^,^^,. (B.3) 



derivative. With the fourth order approximation, 



we found that the deviation from the first law of thermodynamics was, at worst, of order 10~^ for xo,o ^ 3.7 
and even smaller for lower values of Xo.o- 
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Figure 13. Manhattan norm of mismatch vector, plotted as a function of Xo,o- 

This is called the 'mismatch vector' [54]. 

The system of ODEs requires only the values of B and A, and the values and slopes 
of a and X) ^it a single point to completely fix the behavior of the metric and scalar fields 
throughout the spacetime. If the system is reexpressed as a set of six 1st order ODEs, then 
a nonzero mismatch vector means that one or more fields are discontinuous across z^. This 
implies that our initial choice of X produced inconsistent initial conditions at the horizon 
and boundary. The correct choice of X must lead to M = 0. By thinking of M(X) as a 
vector-valued function, the problem becomes that of root finding in six dimensions. We apply 
the Newton- Raphson method (see, e.g., ref. [55]). It works by a generalization of the familiar 
one-dimensional method of tracking tangent lines. For a guess X, compute the Jacobian 
J of partial derivatives of the mismatch vector (Jjj = djMi). Then form the new guess 
Xnew = X — J~^M. Iterate to produce a sequence of X's that converge to the true root. 

We implemented this algorithm in Mathematica 7.0. We let Zmin = d and Zmax = 1-d 
with d = 1/1000. We perform a series expansion of the solution close to the boundary, and 
close to the horizon. Near the horizon the series is evaluated up to and including terms of 
order (z—l)^ for each of the four supergravity fields. Near the boundary the series is evaluated 
up to and including terms of order z^^ In^ z for B, z^^ In^ z for A, z^^ In^ z for a, and z^ In^ z 
for X- The equations of motion were integrated using NDSolve with WorkingPrecision set to 
40 digits, MaxSteps set to oo, and PrecisionGoal and AccuracyGoal each set to 20 digits. 
The matching was performed at = 1/2. For Newton's method we found that a step size of 
1 /lOOOOOO was adequate to compute a forward finite difference approximation to the Jacobian. 
We iterated until the Manhattan norm of the mismatch vector, |M| = Yl^=i l^d) below 
a threshold of 10"'^. In practice, for a sufficiently good starting guess for X, only one or two 
Newton steps were needed to obtain an incredibly small norm. The data in figure 13 shows 
how low our norms became after iterating Newton's method more than a couple times. 

For good numerical performance, it is important to independently integrate inward from 
both the boundary and horizon, since z = and 1 are regular singular points of the ODEs. 
To see this, observe that the field equations for a and x reduce to decoupled Euler differential 
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equations. In all instances, the roots of the indicial equation differ by an integer, and 
this means that one of the linearly independent solutions has a logarithmic singularity. Any 
singular behavior is unwanted near the horizon where the fields should be completely regular. 
Therefore, we seek a solution in which the coefficient of ln(z— 1) vanishes as we approach z = 1 
from below. However, integrating toward this endpoint is a bad idea since the logarithm grows, 
and the solution we do want stays constant. Numerical error has the effect of putting noise 
into the coefficients of the linearly independent solutions, so if the unphysical solution grows 
faster than the good solution, the relative size of numerical error becomes worse the closer 
one approaches to the horizon. It is also unhelpful to integrate toward z = from above since 
both linearly independent solutions vanish and a numerically ill-behaved limiting procedure 
would be needed to isolate the unphysical solution. These practical issues are alleviated by 
starting from the correct solutions at the horizon and boundary, and integrating into the 
bulk where no singular points exist. Although our choice for the "middle of the bulk" was 
arbitrary, we verified insensitivity to its choice a posteriori. 

A crucial question that must be addressed in any numerical computation is: how many 
significant digits in the final result can be trusted? We made extensive use of Mathematical s 
ability to represent arbitrary precision numbers. Testing showed that integrating with 40 digit 
working precision was more than sufficient to guarantee that the mismatch vector's individual 
components did not change to within 10""^. In particular, our results are also sensitive to 
the value of d and the truncation order n of the series. The sense in which d = 1/1000 is 
small must be examined in light of how large n is. There is a simple relation between the 
two: making d smaller is akin to making n bigger since each additional order in a series is 
roughly suppressed by a factor of d compared to the previous order. We kept d = 1/1000 
fixed and decreased n to the next nontrivial order for both series solutions. Solving for the 
roots as before (stopping when a threshold of 10~^ for the norm of the mismatch vector was 
crossed), we found that the new roots were always identical to the original roots in the first 
7 significant digits. 

To generate good initial guesses for the roots, we linearly extrapolated the last two known 
roots along the xo,o-axis. We found 4581 roots spaced at intervals of 0.001 from xo,o = to 
4.58. For the case xo,o = 0; the exact root is X(0) = (— 2, 0, 0, In \/2, 2, 0, 0), which follows 
from the AdS-Schwarzschild solution.'^'' Naive guess-and-check was employed to find the first 
few roots for xo.o > 0- Although Newton's method can diverge rapidly if a poor starting guess 
is made, in practice, we found it to be quite forgiving. Slightly fancier methods involving 
backtracking were used (quite infrequently) to search for roots when the linear extrapolation 

^^Near z = they become a" - {3/z) a' + a = and x" - (3/z) x + (3/^^) X = 0. Near z^l they 

become a" + a'/{z-l) = and x" + x7(-z-l) = 0. 

standard form of the AdS-Schwarzschild metric (in units where L = 1) is ds^ — p^[~f{p) dt^ + dx^] + 
p~^dp^ / f{p), where /(p) = 1 — {ph/pY and pu = ttT. This can be rewritten in the form of eq. (A.l) by 
defining p — pu Vcosh 2r. The horizon is located at r = 0. A further change of variables to z — e"'' puts 
the metric into the form of eq. (A. 3), from which we find that B{z) = (1 — z*)/{l + z'^). A Taylor expansion 
around the origin gives B{z) = 1 — + 0{z^) from which we read off B^.o = —2. Also, Bh = — -B'(l) — 2. 
From eq. (A. 11) we find that Ah = In v^- Lastly, the scalar fields are identically zero in this background. 
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failed to produce a reasonable guess [55] . 



C Fluctuation analysis 



Tensor channel 



For numerical integration it is convenient to consider the system 



xp' + (4A' + B'/B - 3/z)i; + K^e'''^'^(p = 0, 



(C.l) 



Since z = 0, 1 are singular points, we use series solutions to provide initial conditions at ^min 
and Zmax) then integrate inward to a matching point z.^. The regular near- horizon solution 
may be written in the form 



We evaluate the coefficients up to cg. The normalizable near-boundary solution has the form 



where we compute the coefficients up to 68,4. 

Our algorithm for solving the small fluctuation equations is similar to that used for the 
background and discussed in Appendix B. For a given xo,0) the near- horizon amplitude cq is 
arbitrarily set to 1 and a guess is made for the eigenvalue k and the near-boundary amplitude 
60,0 • Next (j) and its derivative V are evaluated with the series at Zmin and ^max, which are 
then fed as initial conditions into NDSolve. Using interpolations of the numerically-generated 
background fields, the system of ODEs is integrated inward to the matching point z*. We used 
machine precision settings for NDSolve. The mismatch vector for ((/>(z,,), is computed 

and Newton's method in two dimensions is iterated (the Jacobian was calculated with a step 
size of 1/1000). Using 5 Newton steps, we obtained Manhattan norms typically several orders 
of magnitude below 10^^. This yields an estimate for the root k in units of 1/L. To find k 
in units of vrT, we multiply the root by 2/{e^^Bh)- 

Vector channel 

The fluctuation equations in the vector channel are 





n 



(Piz) = 60,0 1 + 60,1 In z + ^ ^ 6„,fc(ln z)'' 



(C.3) 



n>l fc=0 



■tp' + {4A' - B'/B - 3/z)^ + K^e-'^'^<p = 0, 



(p' = tp. 



(C.4) 



The regular near-horizon series solution is 




n>l 



(C.5) 
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l + 6o,ilnz + ^z"^6„,fc(ln 



(C.6) 



for which we compute the coefficients up to C5. The normahzable near-boundary series solu- 
tion is 

(fiiz) = 60,0 

n>l fc=0 

for which we compute the coefficients up to eighth order. The numerical algorithm is identical 
to that for the tensor channel. 

Scalar channel 

Derivation of gauge-invariant equations for helicity zero 

We outline the necessary steps to obtain eqs. (4.15). The 3z, 00, 33, zz, and sum of the 
11 and 22 Einstein equations give (in order) 

Hqo + + b'Hoo = -8(3aci5 + x'ciX), (C.7a) 



i/^o + (5^' + 26' + 1/^)^00 + + + H33) + ^^^—Hoo (C.7b) 

8(y„5 + V^x) 



-2A 



Z2 



3^2 



H'^ + (5^' + h' + l/z)H'^^ + A'iH'oQ + H'^) + k2 (^00 + H+) (C.7c) 



8{Va5 + F^x) 



3z2 

H'^, + Hl + H'i, + i2A' + 2b' + l/z)H',o + {2A' + l/z){H'^ + H',,) (C.7d) 



-16(3a[,ia' + x'dX' 



8(y„5 + V^x) 



3^2 

i^- + i6A' + b' + l/z)H'^ + 2^'(i?^o + ^33) + = -l^i^:i|±M) . (c.7e) 



The scalar field equations are 



5" + + 5' + 1/^)5' + ia;i(i/^o + K + ^33) + '^'^S = ^""«^+^"xX ^ ^^^^^ 



r + (4^' + b' + l/z)x' + ix;i(^^o + + H',,) + K^^x = ^""^,t2^""" . (C.Sb) 



,-2A 

2z^ 

Hence, there are seven coupled ODEs in axial gauge. 

In simplifying the above system it is helpful to use the background equations. These 
relations make it possible to express second (or higher order) derivatives of the background 
fields in terms of their ffist derivatives. The ffist step is to recognize that does not 
appear in any of the ODEs, only its derivatives, so one more equation is ffist order in //gg 
and can be considered a constraint. This allows one to drop one of the 2nd order equations, 
say eq. (C.7d). The second step is to eliminate iJgg from the remaining equations. To do 
this, sum eqs. (C.7b) and (C.7e), then solve for i^gg. As the other independent equation, keep 
eq. (C.7b). Now plug in for /fgg in all four equations. The third step is to use the residual 
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gauge invariance that preserves the axial gauge, h^z = 0, to ehminate H^. In axial gauge, 
eq. (4.4), there is a residual gauge invariance parameterized by a vector field r?^ = ri^{z)e~'^^ 
which obeys —'V^rjz — ^ z^Hi-l = 0- As a contravariant vector. 



'^(z) = (Co, Ci, C72, C3 + CzK d/e-2^(^')//, Czz), 



(C.9) 



for arbitrary constants C^.'^^ While this preserves /i^^ = 0, consider what happens to the 
other fluctuation components. In particular, 

-ICzzB^e^'^iA' + B'/B), (C.lOa) 
-ACzzA'e^^, (C.lOb) 



dhoo{z 
Sh+{z 

5h-iz{z 

5a{z 



2e 



2A 



C^k + Czk" dz'e-^^'^^'^/z' -CzzA' 



-Czza^^, 



(C.lOc) 

(C.lOd) 
(C.lOe) 



-CzZXd- 

One can easily verify that the fields defined in eq. (4.14) are gauge invariant. Finally, some 
algebra and judicious rearrangement are needed to obtain eqs. (4.15). 

Near-horizon solution 

The regular near-horizon series solutions are 

Z(z) = Zo + J^Z„(z-l)", 



n>l 



a{z) = ao + ^a„(z-l)' 

n>l 

c{z) = Co + y^Cn(z-l)'' 



(C.lla) 
(C.llb) 
(C.llc) 



n>l 



It is straightforward to determine the coefficients by solving eq. (4.15) order-by-order in z — 1. 
We did this up to fourth order. We find that Zi = —2Zq, Si = — 2ao, and ci = — 2co. Higher 
coefficients are lengthy and we do not write them down. 

Near-boundary solution 

The normalizable near-boundary series solutions are 



Z{z) = z Zpfl + Zo,i In z -^n,fc(ln z) 

L n>l k=0 

n 

a{z) = ao,o + ao,i In z + ^ ^ an,fc(ln z)'', 

n>l k=0 

r " 

co,o + co,i In z + Cn,fc (In z) 

n>l fe=0 



C(Z) = z 



(C.12a) 
(C.12b) 
(C.12c) 



'^If Cz = 0, then the residual gauge freedom simply corresponds to translations in the 4d space (t, x). 



- 39 - 



where we compute the coefficients up to fifth order. Odd coefficients up to and including 05^^ 
and cs^fc are zero. Note that ao,i is not fixed by eq. (4.15) — it is the leading coefficient of 
the non-normalizable solution. We set it to zero by hand. 

Numerical procedure 

The shooting method for the coupled set of ODEs cannot start from 2; = or 1, as 
these are singular points. Therefore, it is necessary to evaluate the series solutions on the 
basis vectors at z = Zmin and Zmax to start the numerical integration. We used settings for 
NDSolve and Newton's method identical to those for the tensor channel. However, in this 
case, Newton's method only needs to be applied in one dimension, making the calculation 
much simpler. 
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